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We scrutinize the real-frequency structure of the self-energy in the superconducting state of the attractive 
Hubbard model within the dynamical mean-field theory. Within the strong-coupling superconducting phase 
which has been understood in terms of the Bose-Einstein condensation in the literature, we find two qualitatively 
different regions crossing over each other. In one region close to zero temperature, the self-energy depends on 
the frequency only weakly at low energy. On the other hand, in the region close to the critical temperature, 
the self-energy shows a pole structure. The latter region becomes more dominant as the interaction becomes 
stronger. We reveal that the self-energy pole in the latter region is generated by a coupling to a hidden/ermion/c 
excitation. The hidden fermion persists in the normal state, where it yields a pseudogap. We compare these 
properties with those of the repulsive Hubbard model relevant for high-temperature cuprate superconductors, 
showing that hidden fermions are a key common ingredient in strongly correlated superconductivity. 

PACS numbers: 67.85.Lm, 74.25-q, 71.10.Fd 


A range of metals show superconductivity below a criti¬ 
cal temperature {Tf}, at which paired electrons acquire a spa¬ 
tial coherence. In conventional superconductors, the pair¬ 
ing mechanism is well described by the Bardeen-Cooper- 
Schrieffer (BCS) theory [1]. However, the BCS theory does 
not straightforwardly apply when the attractive interaction be¬ 
tween electrons is strong. In this case, the electron pairing 
occurs at a temperature higher than Tc and superconductivity 
arises when the electron pairs (regarded as composite bosons) 
go through the Bose-Einstein condensation (BEC) at a lower 
temperature [2]. Such superconductivity in the BEC regime 
has been explored for ultracold fermionic atom systems [3- 
5]. In fact, tightly-bound pairs [6] and an associated pseudo¬ 
gap behavior [7, 8] have been observed for gas in the 
strongly-interacting region. A recent experiment also sug¬ 
gested that the superconductivity in EeSe is in the BCS-BEC 
crossover regime [9]. The preformed pairing has also been 
proposed in the context of cuprate high-Tc superconductors 
[10, 11]. Although it looks unlikely that the preformed pair¬ 
ing solely can explain the whole pseudogap behavior in the 
cuprates [12-15], it may be relevant in a region close to Tc 
around the optimal doping [16, 17]. 

On the theoretical side, the crossover from BCS to BEC 
[5, 18-20] has been intensively studied with continuous [21- 
27] or lattice [28—41] fermion models, both of which give a 
similar phase diagram. The latter is, however, more tractable 
with numerical simulations and allows us to employ nonper- 
turbative methods to study this problem. In particular, the dy¬ 
namical mean-field theory (DMET) [42], which becomes ex¬ 
act in infinite spatial dimensions, is a suitable tool to study 
dynamical properties of the lattice models. In fact, the DMET 
and its extensions have been extensively applied to the attrac¬ 
tive Hubbard model [43-58]. Among the results, for strong 
coupling, a pseudogap in the spectral function has been found 
above Tc [52-56]. Eurther, the pseudogap state was found to 
be separated from a Eermi-liquid metal at weak coupling by a 


first-order pairing transition [44-47]. In light of an electron- 
hole transformation [20], this first-order transition is mapped 
onto a well-known Mott metal-insulator transition in the re¬ 
pulsive half-filled Hubbard model under an external magnetic 
field. Since the Mott gap in the repulsive model is character¬ 
ized by an emergent pole in the self-energy, a self-energy pole 
emerges equivalently in the pseudogap state of the attractive 
model [23, 53, 54], signaling the formation of the preformed 
pairs above Tc. 

It is a nontrivial open problem to determine the fate of the 
self-energy pole below Tc since the sudden appearance of the 
pole in the normal phase appears incompatible with a smooth 
crossover of superconducting properties from weak to strong 
coupling. At weak coupling in the BCS region, the dynamics 
of quasiparticles, described by the low-frequency dependence 
of the self-energies, is well understood. On the other hand, at 
strong coupling in the BEC region, the superconducting tran¬ 
sition is well described by the boson condensation of tightly- 
bound Cooper pairs [20], while the dynamics of the origi¬ 
nal fermions, which is directly relevant to the single-particle 
spectroscopies such as the radiofrequency [59]/photoemission 
spectroscopy, has not been fully addressed. 

In this paper we propose an alternative (“fermionic”) view¬ 
point to the BEC of tightly-bound Cooper pairs. This decrip- 
tion also applies to the intermediate-coupling region where 
the bosonic picture does not hold anymore, bridging the BEC 
limit with the well known BCS region. Our main purpose 
is to reveal that the low-energy quasiparticle dynamics in the 
BEC region is governed by a coupling to a fermionic exci¬ 
tation (hidden fermion) arising from the strong interaction. 
The hidden fermion generates the self-energy pole in the nor¬ 
mal states and persists below Tc, similarly to what has been 
found in the superconducting state of the repulsive Hubbard 
model [60]. As the temperature T lowers further, however, the 
hidden fermion loses the weight and eventually vanishes, in 
marked contrast with the repulsive case. These results set out 



2 



FIG. 1: (Color online) Temperature dependence of (a) the supercon¬ 
ducting order parameter and (b) = E‘“’°(0) - E“°(oo), for 

various values of U. (c) Phase diagram against T and |?7|, with in¬ 
tensity plot of -ZlE™®. 


for a unified view of the superconducting transition of tightly- 
bound Cooper pairs and its relation with the Mott physics. 

We consider the normal and superconducting states of the 
infinite-dimensional Hubbard model on a Bethe lattice 
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which has the advantage to be exactly solvable with the 
DMFT [42]. Here the onsite interaction (7 < 0 is attrac¬ 
tive, t is the hopping integral and /r is the chemical poten¬ 
tial. creates (annihilates) a fermion with spin a 

at site i, and riia = al^aio-- The bare density of states is 
po(w) = with D = 2, which determines the 

energy scale in the paper. In order to avoid possible peculiar¬ 
ities at half filling, we set the fermion density per site to be 
n = 0.8 [61]. 

The DMFT maps this model onto an impurity problem, 
which we solve with a finite-T extension of the exact diag- 
onalization (ED) method [62] elaborated in Supplementary 
Information. The ED solver allows us to study precise real- 
frequency (w) structures of the self-energies, which are of our 
main interest. Eor superconductivity, we calculate the normal 
and anomalous self-energies, E"°‘ and related to the nor¬ 
mal Green’s function by 

G(e, u) = [uj-e + fi- - W{e, (2) 


with 


W{e,uj) 


a; + e - ^ + ' 


( 3 ) 


Since in DMET G depends on particle’s momentum 
only through the hare-dispersion energy e, the den¬ 
sity of states (DOS) is calculated from Eq. (2) as 

Eigure 1(a) plots the superconducting order parameter, 
(flit®*!-)’ against T for various values of U. We define Tc 
as the lowest temperature above which is less than 


FIG. 2: (Color online) (a-d) Various real-frequency properties at 
U = —8 in the normal (T = 0.13 > Tc ~ 0.123) and super¬ 
conducting (T = 0.12, 0.11, 0.06, 0) states. The insets to panels (b) 
and (d) are the enlarged views for T — 0 and 0.06. (e), (f) Change 
of ImE"“ and Re[u; — /r -F E"”(~<^)*] across Tc, where rjo = 0.01 
(see Supplementary Information) is used to sharpen the pole and zero 
features. The an'ows in panel (e) indicate the pole created below Tc 
The thick and thin arrows in panel (f) indicate the pole and zero of 
the plotted quantity, respectively, for T — 0.121, which are created 
in the superconducting state. 


0.01. Thus-determined phase boundary in Eig. 1(c) is consis¬ 
tent with previous works [49-54]. 

Eor \U\ < 4, Tc increases with \U\. In this region, the 
superconductivity follows well the BCS theory; a gap opens 
in the DOS only below T^', for T > E"” shows a Eermi- 

liquid-like behavior, i.e., ImE"” oc (not shown here). 

On the other hand, for | ?71 > 4, in the so-called BEC region, 
Tc decreases with \U\. In this case the normal state displays 
preformed pairs and a gap (called pseudogap) in the DOS even 
above Tc [52-56], as shown in Eig. 2(a) for U = —8 and 
T = 0.13 > Tc ^ 0.123. This gap results from a (thermally- 
smeared) pole in E"“ (because W = 0), which can be seen 
for instance at a; ~ —0.7 for T = 0.13 in Eig. 2(b). The first 
finding in this study is that this pole survives even below Tc, 
symmetrically split with respect to w = 0 by the appearance 
of superconductivity which creates Bogoliubov electron-hole 
branches. These peaks can be for example seen for T = 0.12 
at (jj ~ ±0.9 in Eig. 2(b). The second remarkable finding is 
that E™°(a;) also shows the appearance of these poles at the 
same energies uj = ±0.9, which give considerable dynamical 
structure to ReE™° [Eig. 2(c)] and ImE™° [Eig. 2(d)]. Surpris¬ 
ingly, the weight of these poles diminish as T lowers further 
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[Fig. 2(d)], becoming indiscernible for T < 0.06 (inset). Con¬ 
sequently, ReE^"® is nearly flat at low temperatures [red and 
green curves in the panel (c)]. The value agrees with what is 
expected for a static mean-field decoupling of the interaction 
term, Urii^ni^, in Eq. (1); its anomalous channel is given by 

Despite this drastic evolution of the self-energies (i.e., ap¬ 
pearance of the poles in E™° at and their subsequent dis¬ 
appearance at lower T), the gap in the DOS does not sub¬ 
stantially change with lowering T [Fig. 2(a)] [63]. The singu¬ 
larity generating the spectral gap is in fact transformed from 
the normal self-energy poles for T > Tc to a pole in W, 
i.e. a zero of Re[a; + e — ^ + E"°''(—w)*], for T < 

In order to elucidate this point, we plot in Figs. 2(e) and 
(f), ImE"“ and Re[a; — ^ + E‘'™(—w)*], respectively, in the 
vicinity of Tc [66, 67]. (Compare also the temperature de¬ 
pendences of the peak weights between ImE"”' and ImW in 
Fig. 3.) At T = 0.123 > Tc, ImE"“ shows a single peak at 
u! = ujnor = —0.65, while, as mentioned above, at T < Tc 
another pole in E"™ is created at the particle-hole symmetric 
energy, uj = —Wnoi, as indicated by the arrows in Fig. 2(e). 
The corresponding pole in Re[a; — fJ, + E"°''(—w)*] (arising 
from the pole of E"™(—w)*) is indicated by a thick arrow 
for T = 0.121 in panel (f). This new pole is accompanied 
by a new zero of Re[u; — /r + E"“'(—w)*], as indicated by 
the thin arrow in the same panel. This new zero is the ori¬ 
gin of the strong negative peak in ImW and Im(E"“ -f W) 
[Fig. 3(a)], which yields the spectral gap in the superconduct¬ 
ing state. This is similar to the BCS superconductors in that 
a zero of Re[a; -f e — p -f E‘’“'(—w)*] gives the supercon¬ 
ducting gap. However, as discussed above, the zero in our 
case is tightly connected to the presence of the normal self¬ 
energy poles which are a smooth continuation of those already 
present in the normal state. Thus, the pseudogap and super¬ 
conducting gap are closely related with each other while they 
involve different singularities. 

Another notable finding is that the pole of E"“ and one of 
the poles of W cancel with each other in G [Eq. (2)]. To see 
this property, we plot in Fig. 3 ImE"”, ImW and their sum 
Im(E"°‘' + W) at e = 0 for several temperatures at t/ = —8. 
Surprisingly, in panels (a) and (b), we And that the polelike 
peaks of ImE"” (blue curve) and of ImW (green) cancel out 
and leave no trace in their sum (red) at the peak energy. As T 
lowers from 0.12 to 0.10 [panels (a) to (b)], the polelike peaks 
diminish and eventually below T = 0.06, ImW becomes neg¬ 
ative everywhere so that no cancellation can occur anymore 
[(c) and (d)]. The cancellation is not found of course in the 
BCS region, where no prominent peak in ImE"” is present. 

Turning back to Figs. 2(c) and (d), we see that the appear¬ 
ance of the low-energy pole makes E™°(w = 0) = ReE™°(0) 
depart from the static value, E™°(c») = ReE“°(c») = 
U{ai^aii). We utilize the difference, Z\E™° = E™°(0) — 
E™°(oo), to measure the significance of the self-energy pole 
on low-energy physics. Figure 1(b) plots ZlE™° as a function 
of T for various values of U. In contrast to the absence of any 
appreciable T dependence at weak couplings (|[/| < 4), pro¬ 
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FIG. 3: (Color online) Presence/absence of the cancellation of the 
peaks between ImE"“ and ImW for various temperatures at (7 = 
—8. The insets are the enlarged views of the each panel. Orange 
vertical lines indicate the energy of the peaks which cancel out. 


nounced negative peaks appear at strong couplings (\U\ > 4) 
just below Tc. In Fig. 1(c), we have superposed the intensity 
map of — Z\E™° on the phase diagram. The strong intensity is 
found only on the fringe of the strong-coupling superconduct¬ 
ing phase. This colored “dynamical” area is characterized by 
the presence of the self-energy poles with the canceling prop¬ 
erty, as discussed above. 

As it is shown in Ref. 60, the cancellation of poles suggests 
a presence of a hidden/erm/o«/c excitation (/) hybridizing 
with a quasiparticle (c). For sake of clarity, let us briefly repeat 
a proof of this statement. The hybridization V between c and 
/ would be most simply accommodated in the phenomeno¬ 
logical Hamiltonian, 

i7TCFM = [ec4ca + ^fflfa +V{clU + fUa)] 

(7 

-(TcCtC| + T//tA+h. c.), (4) 


where ec(e/) is the bare one-particle energy of the c(/) 
fermion. For non-zero values of Dc and Df, the system is 
in the superconducting state. By integrating out the / degree 
of freedom in the path integral of the corresponding action, 
we obtain the normal and anomalous self-energies of the c 
fermion as 




e/ 




E“°(w) = D^-V‘^ 


Df 


u;^-e}-Dy 


( 5 ) 

(6) 


These equations show that both self-energies have poles at 
the same energies, w = ±i/ei + D‘f, consistently with 


Figs. 2(b) and (d). By deflning Wc{uj) = 

similarly to Eq. (3), we can calculate the residues of the 
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poles at w = ± y for both and Wc, to obtain 

^ (l ± (l ± 7^) > respectively. 

Therefore, their residues cancel out in their sum. In Supple¬ 
mentary Information, we show that Eqs. (5) and (6) indeed 
fit the DMFT self-energies nearly perfectly at low frequencies 
when the cancellation occurs [Figs. 3(a) and (b)]. On the other 
hand, Fqs. (5) and (6) are not compatible with Figs. 3(c) and 
(d) where no cancellation is observed. The latter is better fitted 
by a fermion-boson model [60]. Thus, the pole cancellation 
gives a stringent test for a fermionic origin of the self-energy 
peaks. 

Interestingly, the model (4) describes well the low-energy 
DMFT results in the normal state, too. In fact, the single pole¬ 
like structure of E"”, seen for T = 0.13 in Fig. 2(b), is con¬ 
sistent with Eq. (5), which yields 

srXw) =- (7) 

uj-ef 

for Dc = Df = 0. 

This observation gives us a clue to identify the hidden 
fermion / in Eq. (4). In the strong-coupling region close to 
the atomic limit, the normal-state Green’s function is well ap¬ 
proximated by 




1 — m 


uj + fj, — U 


u) + ji— Um — 


U‘^m{l — m) 

’ + fi — U{1 — m)_ 


( 8 ) 


with m being the average c-fermion density per spin. Then, 
comparing Eq. (7) with the second line of Eq. (8), we find that 
Gatm is reproduced by the model (4) with V = Uy/ m(l — m) 
and ef = U{1 — m) — fi. 

Thus, the hidden fermion in the normal state, generating 
the pseudogap, is identified with the hidden fermion that ap¬ 
pears by taking the atomic limit of the Hamiltonian (1). It 
is then remarkable that the DMFT self-energy changes con¬ 
tinuously across Tc [Figs. 2(e) and (f)] since the continuity 
suggests that / found below Tc is a smooth continuation of 
the one originating the pseudogap above Tc (where the atomic 
limit is most valid). The local origin of / in turn indicates its 
ubiquitous nature in the strong-coupling region, irrespective 
of dimensionality and lattice structures. 

Thus, we have revealed the presence of the hidden 
fermionic excitation which dominates the quasiparticle dy¬ 
namics in the strong-coupling superconductor just below Tc. 
In order to quantify the area involving the hidden fermion, we 
define Tq as a temperature at which changes its sign 

from positive (for T < Tq) to negative (for T > Tq). Below 
To, the order parameter is almost saturated [Fig. 1(a)] and the 
self-energy becomes nearly frequency independent [Fig. 2(c)]. 

Figure 4(a) and the inset show that the hidden fermion 
indeed takes a major part of the superconducting phase for 
strong couplings. Furthermore, Fig. 4(b) shows a remark¬ 
able scaling, Tq/Tc ~ |( 7 |“ with a ^ —0.36, which indicates 
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FIG. 4: (Color online) (a) Tc and To (defined in the text) plotted 
against |?7|. Inset shows a ratio, (Tc — To)/Tc, of the fermionic 
region to the superconducting region, (b) Scaling of Tq/Tc to \U\. 


that the hidden-fermionic area covers all the superconducting 
phase in the strong-coupling limit. 

Fastly, we compare the above results with those obtained 
in the two-dimensional repulsive Hubbard model. In the lat¬ 
ter, cluster extensions of the DMFT [68, 69] have ascribed the 
origin of the pseudogap and the associated low-energy pole in 
above Tc to “Mottness” emerging in the proximity to the 
Mott insulator [15, 70-79]. In this case too a two-component 
fermion model well describes the pseudogap above Tc as the 
hybridization gap formed with a hidden fermion [80]. Be¬ 
low Tc, particle-hole symmetric poles appear in E"” and E™° 
[81] in a fashion similar to the attractive case, displaying the 
same pole cancellation and a similar mechanism of opening 
the superconducting gap at zero of Re[a; -F e — p -F E'‘™(—w)*] 
[60]. On the other hand, an essential difference is in the en¬ 
ergy scale of the pseudogap, which is of the order of |( 7 | in the 
attractive case while it is less than t in the repulsive case, sug¬ 
gesting different characters of the hidden fermion /. Another 
important difference is in the role of / in the superconductiv¬ 
ity: In the repulsive case / helps to enhance Tc [60, 82], while 
in the attractive case / competes with superconductivity (the 
gap function is slightly suppressed [Fig. S2(a)], even though 
E™° is enhanced by the hidden fermion [Fig. 2(c)]), disap¬ 
pearing at low temperatures where a larger superconducting 
gap is preferred. 

To summarize, we find a hidden/erm/on/c excitation on the 
fringe of the strong-coupling superconducting phase of the at¬ 
tractive Hubbard model. Although the superconducting tran¬ 
sition in the strong-coupling region can be regarded as a bo¬ 
son condensation of tightly-bound Cooper pairs [20], our re¬ 
sults give an alternative fermionic point of view: The quasi¬ 
particles hybridize with a hidden fermionic excitation, which 
originates from a strong-coupling effect, above Tc and down 
to ~ To below Tc. Above Tc, the hidden fermion yields the 
pseudogap in the quasiparticle spectra. It persists in the super¬ 
conducting state down to T ^ Tq, where the atomic physics 
still dominates the quasiparticle dynamics. It disappears be¬ 
low Tq, where the smooth crossover to a static superconduct¬ 
ing state takes place. The hidden fermion region dominates for 
\U\/D —?> oo. Many issues remain open: What is the physical 
entity and the explicit expression of /? Why does E“° lose 
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its frequency dependence below Tq? Why does Tq/T^ follow 
the scaling in Fig. 4(b)? The presence of the hidden fermions 
in two different unconventional superconducting states (for 
U > 0 and U < 0) calls for future search of similar ex¬ 
citations in other unconventional correlated superconductors 
[83-87], 
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SUPPLEMENTARY INFORMATION 
Method 

We solve the DMFT impurity problem with the exact diag- 
onalization method with 8 bath sites for all the presented data 
while we have checked its validity with several calculations 
with 10 and 12 bath sites. 

In order to correctly capture a spectral structure in the 
gapped states, we first arrange the bath sites in a one- 
dimensional-chain form at T = 0 [77]. For calculations 
at T > 0, we rearrange the bath sites so that the im¬ 
purity site couples to every bath sites [78], by fitting the 
hybridization function (obtained by the above chain algo¬ 
rithm) on real-frequency axis. This rearrangement offers 
an initial guess of the bath parameters at low tempera¬ 
tures. We then apply the standard finite-T ED algorithm 
of Ref. 62, with gradually increasing the temperature from 
T = 0. In the algorithm, we optimize the bath parame¬ 
ters at each temperature by minimizing the distance function, 

d = 

to fit the normal/anomalous component of the Weiss func¬ 
tion with that of 8 bath sites Here a;„ = 

(2n + l)7rr is the Matsubara frequency. Note that, while the 
above procedure improves the spectral structure beyond the 
gap energy, the low-energy structure of the self-energy (within 
the spectral gap) is robust against different choices of the ini¬ 
tial guess. 

To display the real-frequency properties, we introduce the 
energy-broadening factor, iri{uj) = ir]Q min[l + 10] to lu. 

We set rjo = 0.05 unless otherwise mentioned. The presented 
results do not essentially depend on the choice of r]{uj). 

Superconducting gap function 


Figure S5 shows the superconducting gap function, which 
is defined through E"" and by [79] 



FIG. 5: The superconducting gap function at [/ = —8 for the same 
temperatures as those used in Figs. 2(a-d) in the main text. 
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FIG. 6: Change of E™° just below Tc ~ 0.123. Inset in panel (a) 
plots the weight of the low-energy peak of ImE™'’ against T. 770 = 
0.01 is used to sharpen the pole and zero features. 


Interestingly, the w-dependent structure of A(a;) is consider¬ 
ably different from that of [Figs. 2 (c) and (d) in the 

main text], suggesting a substantial role of the w dependence 
in E"”(a;). The deviation manifests the departure from the 
BCS and Migdal-Eliashberg theories, which neglect the w de¬ 
pendence of 

In particular, at T = 0.12 where the peak of the hidden 
fermionic excitation is seen at w = ±0.9 in ImE™°, ImA 
does not show any structure at the same energies. Instead, 
ImA shows an intensity at higher energies (|w| > 1.5). 

At T = 0 and 0.06, ImA shows finite amplitudes only in a 
frequency range similar to that of ImE™° while a simple linear 
relation between them (as assumed in the Migdal-Eliashberg 
theory) does not yet hold [see inset to Eig. 2(d) for compar¬ 
ison] as is apparent also from the residual u! dependence in 
E"" [inset to Eig. 2(b)]. 


Continuous development of E^"" around Tc 


nite residue immediately below Tc, where the new zero of 
Re[a; + e — fi + E"°''(—w)*] is created. 


Fitting with the two-component fermion model 


Eigure S7 shows the fitting of the DMET self-energies with 
Eqs. (5) and (6), with Dc = U{c-\-ci) = —0.736 calculated 
by the DMET. We can see that the fitting works nearly per¬ 
fectly for the low-energy peaks at w = ±0.9; Only the visible 
deviation in the presented scale is the weak structures around 
60 
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FIG. 7: Self-energies calculated with the DMFT (red curves) for 
U = —8 and T = 0.12 and with the two-component fermion model 
(TCFM, green curves) for e/ = —0.630, V = 3.39, Df = —0.624, 
and Dc = U{ci-ci) = —0.736. 


Figure S6 shows the continuous development of E^™ 
around Tc, in relation to Figs. 2(e) and (f) in the main text. The 
continuity would be most directly seen in the inset to panel 
(a), which plots the weight of the low-energy peak in ImE“° 
against T. Despite this continuous change of E^™, the singu¬ 
larity responsible to the spectral gap is replaced from a pole 
of E"“ (at uj = Wnoi) to a zero of Re[u; + e — fi + E"“(—w)*] 
(just above Wnor) immediately below Tc as the cancellation 
of the self-energy poles occurs [Fig. 2(f)]. This is possi¬ 
ble because the pole of IE (at w = Wnor) aquires a fi- 


oj = ±2.5 in the DMFT self-energies, which cannot be cap¬ 
tured by the model (4). 
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